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All the recent numerical work on QCD at finite density has used staggered fermions, for well- 
known reasons of economy. In order to get the right number of flavors one has to take the square 
root and the fourth root of the fermion determinant. I will point out a fatal flaw in this procedure, 
one that will invalidate calculations performed at sufficiently large values of /J. unless they are 
performed very close to the continuum limit. The argument in a nutshell is as follows 

The determinant A[U] = Det [D(U) +m] may be written as a product of eigenvalues 17,- of the 
staggered Dirac operator D(U) +m. The eigenvalues are complex, and they move as the gauge 
field U is changed. When Re ;U > the determinant is complex. A reweighting ||] procedure 
needs the fourth root A[£/] 1//4 . One has to define the phase of this root carefully. For instance, if one 
takes the root of each T],- in the cut plane, and defines A 1 / 4 = YlVi > tnen A 1//4 will jump whenever 
any eigenvalue crosses the cut. 

Are these jumps a bad thing? I will show you an alternative prescription for A 1 / 4 that has no 
jumps — and that is obviously correct. So the jumps are indeed wrong, and a source of systematic 
error. Our "ideal" prescription, on the other hand, only works if there is an almost unbroken taste 
symmetry in the low-energy spectrum, which is only true very near the continuum limit. Far from 
this limit, taste symmetry is badly broken and there is no way to keep the jumps from returning, 
along with enormous systematic error. We contend that current calculations are still too far from 
the limit, and that this will be the situation for the foreseeable future. 

I will not touch on the recent progress in demonstrating the validity of the continuum limit for 
staggered fermions at }X = 0. This is the subject of several talks at this conference [Q]. I recommend 
rooting around in the excellent summary by Steve Sharpe [Q]. 

1. Eigenvalues of the staggered Dirac operator 

When /J. = the Dirac operator D[U] is anti-hermitian, so the eigenvalues of D + rn lie on a 
line parallel to the imaginary axis, Re 17 = m (see Fig. [I]). If we are close to the continuum limit, 
the low-lying eigenvalues are near-degenerate because of the approximate taste symmetry; they are 
grouped into quartets corresponding to the four tastes, with an 0(a) splitting. 1 When Re/x > the 
eigenvalues spread out into the complex plane as shown. The density of eigenvalues (per unit area 
in the complex plane) is proportional to the 4- volume V; the splitting of each quartet is still 0(a). 

How wide is the spread from the original Re 77 = m line? The distribution can reach the 
imaginary axis for special configurations U if Re /I > m, the quark mass. 2 Such configurations, 
though, will be rare and atypical until /I is considerably larger. Analysis in the e regime |^] shows 
that the typical eigenvalue distribution will reach the axis (and the origin) when /I > m K /2. This 
prediction is for T = 0; for T > the bound is expected to move upward, as it is connected to the 
phase transition in the phase quenched theory 

2. Fourth roots 

Now let's see where phase jumps come from. We can consider one quartet at a time, and define 
'But see Sec. 5. 

2 Our n is the quark chemical potential, so the baryon chemical potential is flB = 3/1. 
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Figure 1: Sketch of complex eigenvalues 77, of D + m when jj, —0 (left) and when Re fi is somewhat larger 
than m K /2 (right). The dashed line is Re 77 = m, 

tentatively its contribution to A 1 / 4 as a product of roots, each taken with a fixed cut: 

~ 1/4 4 
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The phase of this product will jump by n/2 whenever any r\ crosses the cut (Fig. |||). 



(2.1) 



Figure 2: An eigenvalue quartet as it crosses the cut, which was chosen to lie on the negative real axis. 

These jumps can be avoided by an alternative prescription. Working with one quartet at a 
time, we define the phase of (n^;) 1 ^ 4 to P om t to the center of the quartet, no matter where the 
quartet may be (Fig.^j). Then the phase of the root follows the quartet around without jumps. This 
is nothing other than the smooth replacement of four tastes by one quark per flavor — which is the 
whole idea of taking the fourth root! We call this the "ideal prescription." 

The ideal prescription runs into trouble when applied to certain problematic configurations. If 
there is a quartet that frames the origin (see Fig. [|), there is no clear definition of the phase of the 
root. Whether we keep such configurations or drop them, they will contribute a systematic error to 
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Figure 3: Choosing the phase of the root in the "ideal prescription." Now there is no need for a cut. 



Figure 4: An eigenvalue quartet that frames the origin. 



any result. The relative error will be the probability of finding such a quartet at the origin, which is 
the eigenvalue density times the area occupied by the quartet, or 0(a 2 VA ( ') for quenched configu- 
rations. Reweighting reduces the contribution of these configurations by a factor of 0(aVv A 3 ), so 
the error estimate becomes 0[(a\/V A 3 ) 3 ]. 

The devil is in the factors of volume in these error estimates. They are connected with the fact 
that we must have ayV A 3 <C 1 in order to have well-separated quartets. The volume, however, 
cannot be reduced freely — it is usually fixed by physics requirements, such as m K L > 3 for a ther- 
modynamic "limit." In order to control the error, it is essential to take the limit a — » before taking 

V — » oo. 



3. Other prescriptions 

So much for our ideal prescription. It requires identifying quartets and treating them individ- 
ually. If we don't do this, how bad can things get? 

Let's go back to the fixed-cut prescription, Eq. (2.1). The phase of A 1 / 4 jumps whenever 



an eigenvalue crosses the cut; it is restored to the correct value only when an entire quartet has 
crossed the cut. Thus — if there are clean quartets in the first place — the systematic error comes 
from configurations where a quartet enters an exclusion zone of width ~ a along the negative real 
axis. The systematic error is thus 0(aV A 5 ). (For simplicity we take /I = 0(A); in any case we take 
H > m K /2, so that the eigenvalue density around the origin and to its left is appreciable.) 

If we are not very close to the continuum limit, there will be no quartet structure apparent in 
the eigenvalue distribution. The figure of merit is, again, a\/VA 3 . If there are no quartets then the 
phase flips will be uncorrected and the systematic error incalculable. 
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4. Numbers 

Where are we today? Let's begin with the figure of merit, F = a\/VA 3 . This is the ratio of 
the splitting in a quartet to the distance between quartets. Let's require that it be less than, say, 
1/5. If we fix the 4-volume V =L? /T by setting L > 3m~ l and T ~ T c , and identify the scale as 
A = T c ~ 200 MeV, we arrive at 

a -1 >9GeV, (4.1) 

which should give one pause. This is an order-of-magnitude estimate, which allows either optimism 
or pessimism. The true test in any calculation will be to display the eigenvalue distribution of a 
typical configuration in order to see whether the quartet structure is evident. 

We can judge the state of the art by the largest reweighting calculation done to date [Q], wherein 
N t =4 and T ~ T c , so that cT x ~ 800 MeV, and VA 4 = 27. This gives F ~ 1.3, meaning there can 
be no well-defined quartets to work with. Moreover, the method of Ref . ^ is to write A = ]T 
where the factors ^ are not eigenvalues of D + m. We know of no analysis that predicts a taste 
structure for the ^ and thus no reason why quartets should form in the continuum limit. 

Going a step further, if some calculation does manage to reach a -1 ~ 9 GeV (as it were), then 
F 3 gives an estimate of the systematic error of 0(1%) for the "ideal prescription." Recall that the 
ideal prescription requires identifying all quartets and taking their fourth roots carefully. If, instead, 
one were to use the fixed-cut prescription in this regime in order to avoid identifying quartets, then 
the relative systematic error, calculated as 0(aVA 5 ), would still be ~ 1.8. Bringing this down to, 
say, 10% by decreasing the lattice spacing is a daunting prospect. 

5. a or a 2 ? 

I have assumed that taste breaking in the eigenvalue spectrum is an 0(a) effect, which is the 
weakest assumption consistent with a continuum limit. In discussions following this conference, 
the possibility arose that the splitting of the eigenvalue quartets might be 0(a 2 ) rather than 0(a). I 
add a brief (and inconclusive) discussion of this issue here for completeness. If indeed the splitting 
is 0(a 2 ) then our estimates of systematic errors change in the obvious way; I shall display the 
recalculated estimates below. 

Let us begin with free staggered fermions. At /I = the fermion action may be rewritten in 
terms of block fields, rearranged in the "taste basis," as ^ 

Staste = (2a) 4 £¥(X) [V jU (r M ®/)-aA jU (75«)T M T 5 )+m] y(X). (5.1) 

X 

Here X labels a block of 2 d sites and the fermion field \ff has 16 components, which are decomposed 
into a Dirac index (for / matrices) and a taste index (for z matrices), each running from 1 to 4. The 
first term gives a dispersion relation of naive fermions, replicated for the four tastes. The second 
term, with an explicit a in the coefficient, removes the doublers and splits the taste degeneracy. 
Taking, for instance, = (p,0, 0, 0), with pa -C 1, we find that the first term takes the form ipj\ (g)/ 
while the second is ap 2 y=, (8>TiT5. The two terms don't commute. Diagonalizing the operator 
then gives splittings that are of 0(p 2 a 2 ) relative to the unperturbed energy p. This result, which 
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generalizes easily to any p, shows that the free theory's quartets are split in 0(a 2 ), and not as 
assumed above. 

In the interacting theory, the gauge field insinuates itself into taste-invariant and taste-breaking 
terms alike, and it looks like there is no reason why splittings shouldn't be of the naive 0(a) 
according to the explicit a in the second term in Eq. ( |5.1[ ). (As noted in Ref. [|J], this is not 
inconsistent with the well-known a 2 scaling of all physical observables.) It is known, however [Q], 
that all 0(a) terms in the interacting action can be eliminated by switching to improved fields 
defined by 

W^(l + ^)¥, (5-2) 
where & = 2?(U) connects nearest neighbors and is anti-hermitian. This can be readily general- 



ized to jli / (where 2F is no longer anti-hermitian). This follows [Q, |0|] from the fact that the 
Symanzik effective action for staggered fermions equals the continuum action plus terms of order 
a 2 . What is unclear however, is whether this also implies that eigenvalue splittings are of order a 2 , 
at least for eigenvalues of order A or below (which are the ones at issue here [|T]]). The problem 
is that the usual Symanzik analysis of irrelevant operators applies to the calculation of low-energy 
physical observables, not to the calculation of individual eigenvalues of the Dirac operator. It thus 
remains an open question whether taste splittings in the eigenvalues scale as a or as a 2 . 

If we suppose that the taste splitting is of 0(a 2 ), then the figure of merit becomes F' = 
a 2 WA 4 = (aA)F. Requiring F' < 1/5 gives [cf. Eq. Q] 

a _1 >1.3GeV, (5.3) 

which is no longer astronomical. A reassessment of the parameters used in Ref. pi], moreover, 
brings its figure of merit down to F' ~ 0.3, which may offer a glimmer of hope. 

Two things, however, are unchanged. First, the relative error in the ideal prescription, F' 3 ~ 
1%, applies only if one identifies quartets and treats them as indicated (the Fodor-Katz method, 
as we have mentioned, doesn't do this). Second, if one accepts F' = 1/5 but uses a fixed-cut 
prescription, then the relative error here (a 2 VA 6 ) still turns out to be ~ 1.8, because a is larger than 
before. Even if one starts with cT 1 = 1.3 GeV, this will have to be raised by a factor of y/lS to 
reach an error of 10% (the error, remember, is now ~ a 2 ). 

Whichever error estimate we adopt, it can be reduced by using an improved action. A tree- 
level improvement to the staggered action will reduce taste breaking by a factor of a s , and our 
estimates will change accordingly. We don't see how this will change any orders of magnitude, but 
some improvement can be expected. 



6. Conclusions 

I will end by emphasizing three points. First, the only way to see whether a calculation will 
have its fourth-root error under control is to see whether the eigenvalues calculated for generated 
configurations really fall into well-defined quartets in a sizable region that includes the origin and 
the negative real axis. (More distant eigenvalues can be treated by a fixed-cut prescription.) This 
will render moot any arguments about our estimated figure of merit and about a vs. a 2 . 

Likewise, the whole problem arises only if Re pL is sufficiently large that eigenvalues have 
begun to circle the origin in the complex plane. I have quoted a rough bound Re^i > m n /2 for 
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T = [Q], and noted that this bound should move upwards as T is increased ^J§. Any calculation 
that claims to lie below the bound should prove it by, again, plotting eigenvalue distributions. 

Finally, if the quartet structure is verified, one cannot just go ahead and use a fixed-cut pre- 
scription for all eigenvalues. The error will still turn out large, because there will always be a 
quartet that straddles the cut and gives the wrong phase to the determinant. One must use the 
ideal prescription, proceeding quartet-by-quartet to assign an unambiguous phase to the root of the 
product. 
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